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Casimir forces between conductors at the sub-micron scale cannot be ignored in the design and 
operation of micro-electromechanical (MEM) devices. However, these forces depend non-trivially 
on geometry, and existing formulae and approximations cannot deal with realistic micro-machinery 
components with sharp edges and tips. Here, we employ a novel approach to electromagnetic 
scattering, appropriate to perfect conductors with sharp edges and tips, specifically to wedges and 
cones. The interaction of these objects with a metal plate (and among themselves) is then computed 
systematically by a multiple-scattering series. For the wedge, we obtain analytical expressions 
for the interaction with a plate, as functions of opening angle and tilt, which should provide a 
particularly useful tool for the design of MEMs. Our result for the Casimir interactions between 
conducting cones and plates applies directly to the force on the tip of a scanning tunneling probe; 
the unexpectedly large temperature dependence of the force in these configurations should attract 
immediate experimental interest. 

I. INTRODUCTION 

The inherent appeal of the Casimir force as a macroscopic manifestation of quantum "zero-point" fluctuations 
has inspired many studies over the decades that followed its discovery [T]. Casimir's original result [3j for the force 
between perfectly reflecting mirrors separated by vacuum was quickly extended to include slabs of material with 
specified (frequency-dependent) dielectric response [3]. Precise experimental confirmation, however, had to await the 
advent of high precision scanning probes [3HZ] ■ Recent studies have aimed to reduce or reverse the attractive Casimir 
force in practical applications. In the presence of an intervening fluid, experiments have indeed observed repulsion 
due to quantum [8] or critical thermal j9] fluctuations. Metamaterials, fabricated designs of microcircuitry, have also 
been proposed as candidates for Casimir repulsion across an intervening vacuum |lQj. 

Compared to studies of materials, the treatment of shapes and geometry has remained at a primitive stage. Inter- 
actions between non-planar shapes are typically calculated via the proximity force approximation (PFA) , which sums 
over infinitesimal segments treated as locally parallel plates |llj . This is a serious limitation since the majority of 
experiments measure the force between a sphere and a plate, with precision that is now sufficient to probe deviations 
from PFA in this and other geometries [TH [T3] . Practical applications are likely to explore geometries further removed 
from parallel plates. 

The formalism recently implemented in Refs. [HJ HI] enables systematic computations of electromagnetic Casimir 
forces in terms of a multipole expansion. Using these methods we have been able to compute electromagnetic Casimir 
forces between various combinations of planes, spheres, and circular and parabolic cylinders [T5 20J (see also |21H23| ). 
Both perfect conductors and dielectrics have been studied. However, with the notable exception of the knife-edge 
19J, which is a limit of the parabolic cylinder geometry, systems with sharp edges have not yet been studied 1 . 

In all the above cases, the object corresponds to a surface of constant radial coordinate. In this paper we present 
the first results on the quantum and thermal electromagnetic Casimir forces between generically sharp shapes, such 
as a wedge and a cone, and a conducting plane. We accomplish this by considering surfaces of constant angular 
coordinate. While the conceptual step — radial to angular — is simple, the practical computation of scattering 
properties is nontrivial, necessitating complex, and in places novel, mathematical steps. Furthermore, the inclusion 
of wedges and cones practically exhausts shapes for which the EM scattering amplitude can be treated analytically 2 . 



* Present address: Center for Studies in Physics and Biology, The Rockefeller University, 1230 York Street, New York, NY 10065, USA. 

1 Knife edge geometries have been studied for scalar fields obeying Dirichlet boundary conditions in Refs. 1241 1251 . Wedges and related 
shapes have been studied in isolation 1201 . but computing interactions amongst such objects requires full scattering and conversion 
matrices, which arc first synthesized in this paper. 

2 Morse and Feshbach 27 enumerate six coordinate systems in which the vector Helmholtz equation for EM waves is generically solvable in 
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In light of the technical complexity of the analysis, we first describe our methods qualitatively and summarize our 
results in the next section. The rest of the paper is organized as follows: In Section |III| we present the analysis for 
wedge geometries, where for perfect conductors the electromagnetic field can be parameterized in terms of two scalar 
fields. In Section |IV| we treat the cone, both for scalar fields and for electromagnetism, where the vector nature of 
the field is unavoidable. In each case we not only describe the formalism but also present further figures and results. 
The role of thermal fluctuations is discussed in Section [Vj where we provide the explicit formulae for Casimir forces 
at finite temperature. In Sec. |VI| we argue that the imperfect conductivity of typical metals leads to controllably 
small corrections to the computed Casimir forces. Finally, in an Appendix we discuss the (novel) representation of 
the electromagnetic Green's function that we employ for computations of scattering from a cone. 



II. OVERVIEW AND DISCUSSION 



The conceptual foundations of the scattering approach can be traced back to earlier multiple-scattering formalisms 
[71 |2"9H5T| . but these were not sufficiently transparent to enable practical calculations. The ingredients in our method 
are depicted in Fig. [T] The Casimir energy associated with a specific geometry depends on the way that the objects 
constrain the electromagnetic waves that can bounce back and forth between them. The dependence on the properties 
of the objects is completely encoded in the scattering amplitude or T-matrix for EM waves. The T-matrices are 
indexed in a coordinate-basis suitable to each object, and by polarization. Fig. 1 summarizes the T-matrices for 
perfectly reflecting planar, wedge, and conical geometries. Scattering from a plane mirror gives T platc = ±1 for 
the two polarizations, irrespective of the wavevector k. Cylindrical (to, k z ) and spherical (t, to) quantum numbers 



label appropriate bases for the cylinder and sphere. As detailed in Sections |III| and IV imaginary angular momenta, 
labeled by fi for the wedge and A for the cone, enumerate the possible scattering waves. The cone's scattering waves 
are also labeled by the real integer to, corresponding to the z-component of angular momentum. The corresponding 
T-matrices depend on the opening angle 9q. The result for the wedge (top-right in Fig. [I]) is independent of the axial 
wavevector k z . In addition to the usual polarizations, for the cone we needed to introduce an extra "ghost" field 
(labelled Gh) and the corresponding T-matrix (top-left in Fig. [IJ . 




FIG. 1: Ingredients in the scattering theory approach to electromagnetic Casimir forces. See the text for further discussion. 



this way: planar, cylindrical (comprising circular, elliptic, and parabolic), spherical, and conical. Surfaces on which one such coordinate 
is constant are candidates for exact analysis. In addition to the plane and sphere discussed above, the circular 28 (bottom-right of 
Fig. Q and parabolic [19. cylinder have already been studied. These shapes are generically smooth, with the extreme limit of the 
parabolic cylinder (a "knife-edge") a notable exception. A survey of the remaining coordinate systems [27| only leads to shapes such as 
cylinders and cones with elliptic cross-section, which are generically similar to their circular counterparts. 
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We will also need the matrix U that captures the appropriate translations and rotations between the scattering 
bases for each object. This matrix encodes the objects' relative positions and orientations. The W-matrices needed 
for our calculations are presented in Sections |III| and |IV| The expression for the Casimir interaction energy, 



he r°° he r°° 

£ = — dK\xM\-m = ~l dn 
2tt J 11 2irJ 



trTV + ItrAf 2 + 



(1) 



involves integration over the imaginary wave number k, an implicit argument of the above matrices, which are 
combined into TV = UT oh ' iCCt U ] <T? utc where we considered one of the objects to be an infinite plane. An expansion 
of tr ln[l — TV] in powers of TV corresponds to multiple scatterings of quantum fluctuations of the EM field between 
the two objects; the trace operation sums over appropriate bases (plane waves for example). This procedure can 
be generalized to multiple objects, with the material properties and shape of each body encoded in its T-matrix. 
Analytical results are restricted to objects for which EM scattering can be solved exactly in a multipole expansion, a 
familiar problem of mathematical physics with classic applications to radar and optics. 

For a perfectly reflecting wedge, translation symmetry makes it possible to decompose the EM field into two 
scalar components: an E-polarization field that vanishes on its surface (Dirichlet boundary condition), and an M- 
polarization field that has vanishing normal derivative (Neumann). In the cylindrical coordinate system (r,faz), 
a wedge has surfaces of constant <j> — ±#o- Whereas for describing scattering from cylinders a natural basis is 
e ±im ^ Hm\i-\/ k 2 + k 2 z r)e %k * z with Bessel-if ^ functions indexed by m = 0, 1,2, for a wedge we must choose 
e ±^07jW (j^/ K 2 _|_ k 2 r)e ikzZ with real \x > 0, corresponding to imaginary angular momenta that are no longer quantized 
(see Section III I. The T-matrices, diagonal in /i, take the simple forms indicated in Fig.[l] Dimensional analysis indi- 
cates that the interaction energy of a wedge of edge length L at a separation d from a plane is £ = — (hcL/d 2 )f(9o, fa), 
where f(9o,fa) is a dimensionless function of the opening angle 9q, and inclination fa to the plane. This geometry, 
and the corresponding function J(9q, <f>o), are plotted in the middle panel of Fig. [5] 

The limit 9q — > corresponds to a "knife-edge," which was previously studied as a limiting form of a parabolic 
cylinder |19) . The matrix TV for the wedge simplifies in this limit, enabling exact calculation of the first few terms in 
the expansion of tr ln[l — TV] in Eq. ([IJ, 
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The first square brackets corresponding to trTV (depicted by an orange line in Fig. [2]) and the second to trTV 2 /2; 
their sum (depicted by a red line) is in remarkable agreement with the full result (blue surface). As fa — > 7r/2, the 
knife-edge becomes parallel to the plane and the interaction energy diverges as it becomes proportional to the area 
rather than L. For parallel plates we know that the terms in the multiple-scattering series tr [TV™] /n are proportional 
to 1/n 4 pQ. Numerically, we find that the convergence is more rapid for fa < ir/2, and that the first two terms in 
Eq. ([2]) are accurate to within 1%. Including more than three terms in the series will not modify the curve at the 
level of accuracy for this figure. Casimir's calculation for parallel plates gives an exact result at fa = tt/2, marked 
with an x . 

Fig. 2 displays some of the more interesting aspects of the wedge-plane geometry. In the middle panel the energy, 
rescaled by an overall factor of hcL/d 2 , and multiplied by cos (6> + fa), is plotted versus #0 and fa. The factor of 
cos (6*o + fa) is introduced to remove the divergence as one face of the wedge becomes parallel to the plane and the 
energy becomes proportional to the area rather than just the length of the wedge. The blue surface is obtained by 
numerical evaluation of the first three terms in the multiple-scattering series of Eq. ([!]) , with TV constructed in terms 
of the plate and wedge T-matrices given in Fig. [I] The front curve (9q = 0) corresponds to a knife-edge as previously 
mentioned. The top panel depicts the "butterfly" configuration where the wedge is aligned symmetrically with respect 
to the normal to the plate. As the wings open up to a full plane at #0 = tt/2, the result again approaches the classic 
parallel plate result, again marked by a x. Finally, and perhaps most interestingly from a qualitative point of view, 
the bottom panel depicts the case where one wing is fixed at 7r/4, and the other opens up by i/j = 29q. This case 
displays the sensitivity of the Casimir energy to the back side of the wedge, which is hidden from the plate. In the 
proximity force approximation the energy is independent of the orientation of the back side of the wedge, and thus 
the PFA result (solid line) is constant until the back surface becomes visible to the plane. The correct result (dotted 
line) varies continuously with the opening angle and differs from the proximity force estimate by nearly a factor of 
two, showing that the effects responsible for the Casimir energy are more subtle than can be captured by the PFA. 

Computations for a cone — the surface of constant 9 in spherical coordinates (r, 9, fa) — require a similar passage 
from spherical waves labeled by (I, m) to counterparts of imaginary angular momentum, I — >■ iX — 1/2. In this case 
A > is real, while m remains quantized to integer values. However, unlike the wedge case, the EM field can no longer 



be separated into two scalar parts; the more complex representation we report in Section IV involves an additional 
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FIG. 2: The Casimir interaction energy of a wedge at a distance d above a plane, as a function of its semi-opening angle 9o and 
tilt <j>o. The rescaled energy as a function of #o and 4>o is shown in the middle panel. The symmetric case, 4>o = 0, is displayed 
in the top panel and the interesting case where the back side of the wedge is "hidden" from the plant is shown in the bottom 
panel. See the text for further discussion. 
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FIG. 3: Casimir interaction of a cone of semi-opening angle 9o a distance d above a plane, scaled by the dimensional factor of 
Fic/d. The left panel corresponds to a vertical orientation, with the energy multiplied by cos 2 9o to remove the divergence as the 
energy becomes proportional to the area for 9o = n/2. The right panel shows the force, F, suitably scaled, for a tilted, sharp 
cone (#o — ¥ 0, evocative of an AFM tip) as a function of tilt angle /3 and temperatures T=300, 80, and OK (top to bottom), at 
a separation of 1/xm. See the text for further discussion. 



field, similar to the "ghost" fields that appear in some quantum field theories. To our knowledge, this is a novel 
representation of EM scattering, which should also be of use for describing reflection of ordinary EM waves from 
cones. Dimensional analysis indicates that for a cone poised vertically at a distance d from a plane, the interaction 
energy scales as (hc/d) times a function of its opening angle 9q. This arrangement and the resulting interaction energy 
are depicted in Fig. [3] (left panel), with the energy scaled by cos 2 9 to remove the divergence as the cone opens up to 
a plane for 8q — > ir/2. The PFA approximation [TT], (depicted by the dashed line) becomes exact in this limit, but it 
is progressively worse as #0 decreases from n/2. In particular, it predicts that the energy vanishes linearly as #0 — > 0, 
while in fact it vanishes as 

_ ftc ln4- 1 1 

d 16tt \\n d f I ' 1 ' 

where the logarithmic divergence is characteristic of the remnant line in this limit [55 . The EM results shown as the 
bold (blue) curve in Fig. [3] are obtained by including two terms in the series of Eq. ([I ) , with Af constructed from the 
plate and cone T-matrices in Fig. [T] including the additional "ghost" field. We have also included the corresponding 
curves for scalar fields subject to Dirichlet and Neumann boundary conditions which are depicted as fine (red) curves 
where the top (bottom) one corresponds to Dircihlet (Neumann) boundary condition. The limit of #0 — ^ is shown in 
the left panel of Fig. 3 as an orange dashed line. As the sharp tip is tilted by an angle /3, the prefactor (In 4 — 1)/I67r 



is replaced by g(/3)/ cos/3, where g(f3) can be computed from integrals of trigonometric functions (see Section IV I. Wc 
plot this quantity in the right panel of Fig. [3j 

For practical applications, the above results have to be corrected for imperfect conductivity and finite temperature. 
The latter correction is easily incorporated by replacing the integral in Eq. (1) with a sum over frequencies k„ = 
(2nkQl '/hc)n. The (analytical) result for the sharp cone is reported in Section|Vj and plotted for d = I/im in Fig. [3] 
Interestingly, the room temperature force is more than 100% higher than T — 0. In contrast, the corresponding 
increase for parallel plates at d = 1/xm is only about 0.1%. This enhanced role of thermal corrections for specific 
geometries has been noted before [32] , and appears essential to the design of MEM devices. Fortuitously, the increased 
importance of thermal corrections diminishes the effects of imperfect conductivity. While we do not compute forces 



for a general frequency-dependent dielectric response e(w), we argue in Section VI that for typical metals (e.g. Au or 
Al), even at zero temperature and for sharp cones (the cases where these corrections are largest), the corrections due 
to imperfect conductivity are at most around 5% for d = 0.2^im. 

Thus for separations 0.2/^m< d< 10/im, relevant to MEM devices and experiments, the perfect conductor results, 
with the important finite temperature corrections, should suffice. For example, let us consider the tip of an atomic 
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force microscope (AFM): At separations, d « 0.2/zm, where the tip may be well approximated as a metal cone, our 
results predict a force that is a fraction of a pico-Newton. Such forces are at the limit of current sensitivities [33], and 
will likely become accessible with future improvements. Current experiments are performed on spheres of relatively 
large radius R, where the force is greater by a factor of (R/d) (a typical R is 100//m). A rounded wedge with radius 
of curvature R falls in an intermediate range, with forces larger [19| by y 1 Rjd than the sharp case. The force can also 
be enhanced by using arrays of cones or wedges, at the cost of the difficulty of maintaining their alignment. 



III. WEDGE 



The Casimir energy and energy density for a wedge in isolation have been considered previously in many contexts 
[3"4Tl4"5] . but these works do not address the interaction energy that leads to the Casimir force. We first formulate the 
scattering theory for a single wedge and include interactions in the subsequent parts. 

For a perfect conductor that is translationally invariant in one direction, it is possible to decompose the EM field into 
two scalar fields that obey the Helmholtz equation with Dirichlet and Neumann boundary conditions, respectively, 
on the conducting surfaces. As the scattering depends only trivially on the translationally-invariant direction, we 
begin by studying the scalar field in two dimensions. In plane polar coordinates the solutions are labeled by the 
component of angular momentum out of the plane, m, and wave number, k. The Helmholtz equation is second-order, 
and therefore has two independent solutions, which we take to be the regular Bessel function J m (kr), which is finite 
at r — 0, and the outgoing Hankel function Hm\kr), which is irregular at r = 0, and obeys outgoing wave boundary 




FIG. 4: The configuration of a tilted wedge opposite an infinite plate. 



conditions for r — > oo. Our primary tool is the free Green's function in polar coordinates for imaginary wave number 
k = in. For applications to scattering theory, a useful representation of this Green's function is in terms of regular 
and outgoing waves [27] . 

oo 

G (r,4>-y,^,iK) = 7 - V e im ^-^')j m («r < )ir ro (Kr > ), (4) 

Z7T * — ' 



in which r< (r>) is the smaller (larger) of r and r'. For imaginary wavenumber, the regular solution becomes a 
modified Bessel function of the first kind J m («r), while the outgoing solution becomes a modified Bessel function of 
the third kind K m {nr). The former diverges for r — > oo, while the latter diverges at r = 0, but Eq. Q avoids these 
pathologies by selecting the regular solution for the smaller value of r and the outgoing solution for the larger value of 
r. Along with the Green's function, we also require the expansion of a plane wave in terms of our scattering solutions, 

oo 

e ifcrcos0 = i m e i,n *J m (fcr) , (5) 



where </> is the angle between k and r. The above representation is suitable for scattering from an ordinary cylinder, 
for which the boundary condition is imposed along a fixed value of r. For the wedge, however, the boundary is instead 
defined by a constant value of <j), so we want the complementary representation in which the discontinuity in the 
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FIG. 5: Analytic continuation to the imaginary axis in 
the complex angular momentum plane for the wedge. 



representation of the Green's function is implemented through 4> rather than r. To create this representation, we let 
the polar angle <f> be defined from — 7r to ir, with +7r and — ir identified. The points at <f> = and cf> — ±7r then serve 
as the analogs of r = and r = oo in defining regular and outgoing solutions. The symmetry axis of the wedge is the 
half-line = and we let 6q be its half-opening angle, see Fig. [4| 

We first carry out this transformation for the expansion of a plane wave. From Eq. ([5| with kr — ix, we obtain 



g — X COS 



= J2 2/ m (z)(-l) m cos(m</>), (6) 



m=0 



where the prime on the sum indicates that the first term is weighted with a factor of 1/2. We can represent the sum 
over m as a contour integral along C of an integrand with poles at non-negative integer values of m. Then the sum 
in Eq. ^ becomes 

e-°-*= [ gL.JL_2I v (x)caaM), (7) 
J c Am sin vn 

where the integration contour is shown in Fig. [5] and the factor of 1/ sin^7r introduces poles with the correct residues. 
The functions in the integrand are analytic, so we can deform the contour to an integral along the imaginary axis plus 
a semi-circle at infinity, which does not contribute to the integral. So we are left with the integral on the imaginary 
axis (C), see Fig. [5] Using 

Ku(x) = - . , 8) 

2 sin v-k 

and the reflection symmetry of Bessel K functions, we find 

2 f°° 

" os ^- / dXK a (x)cosh(X<f>). (9) 
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This integral is convergent only for \<j)\ < 7r/2, since Ki\(x) asymptomatically goes as e A7r ^ 2 for large A. 
More generally we have, for arbitrary complex angle a, 

2 r°° 

e -xcosW-a) = _ / dXK lX {x)cosh(X(cj)-a)), (10) 
n Jo 

which is convergent if |Re(0 — a)| < tt/2. We will consider <j> and a as the angles of r and k respectively. 
The analogous representation of the Green's function can be obtained in a similar way |46j , 

1 f°° 

G {r,cj>;r', ( t>',iK) = ^ dXK iX (Kr)K iX { K r') cosh(A(7r - \<f> - cj>'\)) . (11) 
Jo 

One can also obtain this result using the orthogonality condition on Bessel K functions with respect to their argument, 



dXXsinh(Xn)Kix(Kr)K a (Kr') = r5(r - r') , (12) 



which forms the basis for the Kontorovich-Lebedev transform [47] . 

For scattering theory applications, it is advantageous to cast the Green's function into the same bilinear form as 
Eq. @, 



Go(r,4>;r',(f>',iK) 



Kix(Kr') cosh (A(tt - \<f>>\)) K iX (Kr) cosh(A0<) 



+ sgn^^Kixinr') sinh (A(?r - |0>|)) K iX {nr) sinh(A(^<) 



(13) 



where 0< (</>>) is the angle with smaller (larger) absolute value (|</>|, |</>'|). From this expression we can now read 
off the "regular" and "outgoing" solutions, which are defined as the function which are well-behaved at (j> — and 
<p = ±7r respectively. From Eq. (13 1, we have 



$ r °s< c (r» =i^ A Mcosh(A<£), 
$out,c (/) ^ = XiA ( K /) cosh (A(7r - \4/\)), 

$ ro ^ s (r» =^ A (/ t r)sinh(A0), 
$out )S(/) ^ = XiA ( Kr ') sima (A(7r - |0'|))sgn(</>') , (14) 

where the superscripts c and s stand for symmetric and antisymmetric wavefunctions in (f> respectively. Note that 
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FIG. 6: Plane polar coordinates appropriate for the wedge. The discontinuity at <f> = ±7r constrains the definition of regular 
and outgoing functions. 

the regular functions are irregular at <f> = ±tt. Also the outgoing functions are not regular at <f> = 0. We also define 
the normalization coefficients 



n c/s _ 1 



(15) 



which becomes important upon changing to a different basis, see Eq. (23) 

it ions on 1 
zh we can 

cosh(A6>o 



Dirichlet and Neumann boundary conditions on the wedge are satisfied by an appropriate linear combination of the 
regular and outgoing functions, from which we can read off the T-matrix elements, 



rpC V " / rpS 

DA ~ cosh(A(^-0o))' DX 



sinh(A6'o) 



sinh(A(7r- 6» )) 
cosh(A6> ) 



c sinh(A0 o ) s 

Tnx = sinh(A(7r - 9 j) ' TnX = cosh(A(^ - 0„)) ' 



(16) 



where 9q is the half-opening angle of the wedge. Here D and N stand for Dirichlet and Neumann boundary conditions 
respectively 3 . 



3 The Dirichlet and Neumann boundary conditions in the case of scalar field correspond to electric (E) and magnetic (M) modes for 
electromagnetism, respectively. (D, N) and (c, s) are referred to as (E, M) and (+, — ) in Fig 1. 
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Now let us return to a wedge with half-opening-angle 9q opposite an infinite plate. We take <fi = to be the 
symmetry axis of the wedge. To begin with, we consider the case where this axis is perpendicular to the plane. Then 
the y' axis, defined as the perpendicular axis to the plane, will be parallel to y, the axis along the = line. The 
wavevector, k — (k x >, k y >), satisfies the on-shell condition, hence 



k x , + k y , — k 



(17) 



Of course, k x i and k y i cannot both be real. In fact, k x > , the component of the wavevector parallel to plane, should be 
real while k y / = i\J k 2 + k 2 x , is imaginary. So the angle a as defined by 

k = (k x > ,k v >) — (k x i, i(n 2 + k 2 ,) 1 ^ 2 ) — (zKsina, in cos a) (18) 

is imaginary, and whose range is (— ioo,+ioo). Note that Eq. |l0| ) is convergent for all values of e (— 7r/2,7r/2) (of 
which the wedge is a subset) and can be more conveniently written as 



gifc-r ^ — kt cos(0— a) 



d\ (Ki\{nr) cosh(A^) cosh(Aa) + Ki\(nir) sinh(A0) sinh(Aa)) . 



(19) 



The conversion matrix elements |14| between the wedge and the plane scattering waves are then given by 4 



D 



A.fe x = - cosh(Aa) 

7T 



D\, kx = - sinh(Aa) 



If the plane is tilted by an angle </>0i we have instead 



2 2 
D \,k* = - cosh(A(a - D s x k = - sinh(A(a - <j> Q )) 

VT 7T 



(20) 



(21) 



where the range of a is still given by (— ioo, +ioo). 

So far, we considered the two-dimensional problem and ignored the axis z along which the geometry is translationally 
invariant. Now we rename the two dimensional wave number in the argument of the Bessel functions, which we 
previously called k, to be p, and use k to denote the (imaginary) wave number (u = inc) of the three-dimensional 
problem. So, p — n 2 + k 2 with k z being the wave number in the transverse direction. Because of the translational 
symmetry, we can make a change of variable to eliminate k and recast the three-dimensional problem in a two- 
dimensional basis; the Casimir energy takes the form [HI [15] 



£ = 



HcLz 
2ir 2 



/>oo 

/ dpptrln(l- T1U12T2U21) ■ 
Jo 



(22) 



Here L z is the length of the objects in the transverse direction, U12 and U21 denote the translation matrices between 
objects and T\ and T2 are the objects' individual T-matrices. Note that the "tr" in the last equation is integrating 
over all quantum numbers in the two dimensions (the third direction being included in p) . Using the T-matrix of the 
wedge from Eq. (16) and the conversion matrix of Eq. (21 1, we can express T-matrix of the wedge in the plane- wave 
basis as 



A=c,s 



(23) 



where Ch 



is the normalization coefficient in planar- wave basis as defined in Ref. [14 . We define a = ia 



so that a is real, and convert the k x index in Eq. (23 1 to an a index. We then take a as the outer index in Eq. (22 1, 



with the A index contracted in the matrix multiplication, so that the energy takes the form 

dpptr\n(2 aia > -J\f a>a >) , 



£ 



He L z 
2^T 



(24) 



4 The x-index on k x is unprimed because the conversion matrices are defined with respect to the (x, z) axes of the wedge. 
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where 



KfD — pd(cosh a+cosh a') 



'cosh(A(ia + fa)) cosh(A(ia' — 0o)) cosh(A#o) 



for the Dirichlet case, and 



Af 



N 



-pd(cosh a+cosh a') 



dX 



sinh(A(za - 


cosh(A(7r — 8 )) 
h fa)) sinh(A(ia' — < 


/» o ))sinh(A0o) 




sinh(A(7r - 6» )) 




sinh(A(za - 


- fa)) sinh(A(ia' — (, 


i o ))cosh(A0 o ) 


cosh(A(ia - 


cosh(A(7r — 6*0)) 
f fa)) cosh(A(«a' — 


fa)) sinh(A0 o ) 


sinh(A(7r - 6» )) 



(25) 



(26) 



for the Neumann case, in which d is the distance between the plane and the edge of the wedge and the trace is defined 
as tr f(a, a') = daf(a, a). 

The integral over A can be evaluated exactly, but we first consider a special case. In the limit that the wedge 
becomes a knife edge, i.e., 9 = 0, the second term in the A/"-matrix vanishes and we find 



ktd/n 



3=0 



1 e - P d(co S ha+coshc/) f ±sec i^ a + a >\/ 2 ) + sec(i(a - a')/2 + fa)) 



(27) 



We note that we can also compute the tr In with A rather than a as the outer index. The A/"-matrix in this basis is 
given by 



M D/N 



A, A' 



1 



1=0 



7T cosh(A7r 



(±Jf i(A+A /) (2pd) cosh((A - A')0o) + tf<(A-A') (2pd) cosh((A + A')0o)) , (28) 



and the trace is now defined by tr /(A, A') = J °° dXf(X, A). 

In either description, the indices of the matrix Af are continuous, not discrete. In matrix multiplication, a continuous 
index is defined simply by replacing the sums by integrals. Computing the tr In (or, equivalently, In det) of a continuous 
matrix can also be done by discretization of the continuous index, but for our purposes another approach will be more 
efficient: We will expand trln(l — Af) as a power series in Af and show that we can achieve a sufficient precision by 
keeping only the first few terms. 



The results for a tilted knife edge opposite an infinite plate are plotted in Fig. 7(a). This plot is in agreement with 
the results obtained for the knife edge as the extreme limit of a parabolic cylinder [19 . They have been computed by 
considering the first three terms in the expansion of the tr In in Eq. ( 24 ) . The rapid convergence in powers of Af is 



shown in Fig. 7(b), using the Dirichlet case as an example. For a knife edge (#0 = 0), we can exactly compute the first 
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7(a) Dirchlet, Neumann, and EM boundary conditions. 
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7(b) Convergence of multiple-reflection expansion for 
Dirichlet boundary condition. 



FIG. 7: Tilted knife edge. 



two terms in the expansion of trln formula (i.e., the expansion in multiple reflections). First, we will report results 
separately for the Dirichlet and Neumann cases (denoted by D and N respectively). The electromagnetic result is 
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the sum of the two. We obtain results that are accurate to about one percent by keeping the first two terms in the 
expansion. Using dimensional analysis, the energy is proportional to hcL z /d 2 , so we calculate the scaled energy 



knife 
cedge 
°D/N 



±- 



HcLJd 2 16tt 2 V 8 2 



± 



sec <pQ 



32tt 3 



sec 2 (</> /2) 



1 

12 



16 



(esc (fro sec0o(2^o — sin2</> )) 



(29) 



where the first line is obtained by evaluating trAf and the expression in the second line by evaluating ^trA^ 2 . The 
dots indicate corrections from higher order reflections. It is interesting to study these results in a few different limits. 
When the knife edge is perpendicular to the infinite plate ((fro = 0), the previous equation reduces to 



£ 



knife 
edge 

D/N 



hcL z /d 2 



= ±- 



17 



1 



"2048tt 32tt 2 



192tt 3 



(30) 



a result that heretofore was known only numerically 19 . In contrast, the proximity force approximation (PFA) 
vanishes for this configuration. On the other hand, in the limit where the knife edge is almost parallel to the infinite 
plate ((fro = ir/2 — e), we have 



knife 
r odge 

'D/N 



hcL z /d 2 



32tt 2 



1 

16 



± 



1 



10247T 1927T 3 



0(e)+- 



(31) 



6 o = (ir/2-<0 



The first term diverges as e — > because it gives the contribution proportional to the area that arises as the plates 
become parallel. The second term then represents the leading edge correction which is in reasonable agreement with 
Ref. [TH] for Dirichlet and Neumann boundary conditions. However, the electromagnetic result that sums the two 
varies very slowly (see Fig. 7(a)) and thus the edge term (which is sensitive to the relative accuracy of nearby points) 
deviates from the result in Ref. [19;. We expect convergence to the exact result by including higher-order terms in 
multiple reflections. 

As noted above, the expansion of the logarithm in M corresponds to a multiple reflection expansion. This identifi- 
cation forms the basis for the optical approximation to the Casimir energy [48]. For parallel plates, the terms in this 
series fall in magnitude like 1/n 4 , where n is the n umb er of reflections (back and forth) between the objects. We can 
check this result explicitly for the area term. Eq. (31) displays the first two terms, 1 + in the 1/n 4 series which 
sums to C(4) = 7r 4 /90. The contribution of the first two terms captures more than 98% of the exact result. For a 
wedge (or a knife edge) away from (fr = n/2, we find numerically that higher reflections fall off even more rapidly 
than 1/n 4 . For as many as six reflections, for a given (fro < tt/2, the fall-off is a good fit to l/n 4+5 ^°) where 6 ((fro) is 
a positive function of the angle. 

Finally, the electromagnetic Casimir energy including the first two reflections, is 



knife 
c edge 
'EM 



1 



1 



■ sec (f>o 



+ esc 3 (fro sec (fro (20o — sin20 o ) I + 



HcL z /d 2 16tt 2 ru 256tt 3 V3 
When the opening angle of the wedge is nonzero, the A/"-matrix takes a more complicated form, 



(32) 



■A/" ^ a, a' 



J ,-Mcosha+cosha') x ± aec f '< '" f " ! 



sec 
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2(tt - 9 ) 



sec 
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2(tt - 9o) 
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2(tt~9o) 



± sec 



,'ir(ir + ia + ia')\ ( n(ir + ia — ia' + 2(fr ) 

± COt r — COt ; r 

2(n-6o) J \ 2(n-8o) 

' tt(tt — ia + ia' — 2/fro) \ / tt(tt — ia — ia') 
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2(tt - 9o) 



2(tt-0 o ) 



(33) 
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FIG. 8: Vertical wedge. 



or, in the A basis, 



(34) 



The energies for a wedge of opening angle #0 positioned vertically above the plane, i.e., with tilt angle 0n = 0, 
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FIG. 9: Casimir energy of a wedge with one face at a fixed angle (j>o + So as a function of the full opening angle if) — 28q. The 
PFA prediction, in dashed lines, is compared to the exact calculation, indicated by solid circles, where 3 multiple reflections 
have been taken into account. The PFA does not change until the top face becomes visible to the plane (ip = 2(9o > 4>o + #o)- 
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are plotted in Fig. 8(a), including three terms of the reflection expansion. The convergence in multiple reflections is 
shown for the Dirichlet boundary condition in Fig. |8(b). 

As in the case of the knife edge, we can find analytical results for the wedge (with non-zero opening angle). Only 
including the first term in multiple reflections, the electrodynamic Casimir energy is given by 



p wedge 



1 



hcL z /d 2 16tt 2 2(1 - 6» /vr) 



tan 



7r(7T - 2</> ) 
4(tt - 9 Q ) 



tan 



7T(7T + 20Q ) 
4(7T - O ) 



(35) 



where the dots indicate higher reflections. For Dirichlet and Neumann boundary conditions, the expression for the 
energy is more complicated because the terms of opposite sign in Eq. (331 do not cancel. These terms do not 



contribute to the torque, however, which in both cases is simply one half of the derivative of Eq. (35) with respect to 
4>o- Hence, to this order, the torque is the same for Dirichlet and Neumann boundary conditions. 

The geometry of the wedge provides an interesting example to examine the PFA prediction. Within PFA, the energy 
is computed by integrating over the surfaces facing each other, so that the energy remains constant until the back 
surface of the wedge becomes visible to the plane. Our result of Eq. (35), on the other hand, depends smoothly on the 

0o) 



angles. To demonstrate the failure of the PFA, we consider a wedge with one fa ce fixed, at an angle n/2 — (< 
with respect to the plane, while the other face opens up, as shown in Fig. 9(d). | The energy as a function of the (full) 
opening angle ip = 29q is shown in Fig. 9(a)-[ c). 



IV. CONE 



A. Scalar field 



We now apply the same techniques to the case of a conical perfect conductor opposite a conducting plate. We start 
with spherical coordinates with the origin at the tip of the cone and the z-axis aligned to the cone's symmetry axis. 
The expansion of a plane wave in terms of spherical Bessel functions with imaginary wavenumber reads 



-KT COS @ 



]T(2n+l)(-l)"P n (cos©K( 



KT) 



(36) 



n=0 



where i n is the modified spherical Bessel function. By turning the summation into a contour integration with 
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FIG. 10: The cone coordinates. The limits of the angle 6 are essential in defining the regular and outgoing functions. 
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appropriate poles on the real axis (see Fig. 11 ), we find 

dX Atanh(A7r)fc, iA _ 1/2 (Kr) J R tA _ 1/2 (cose) 

/ dX Atanh(A7r)(-l) m /c 2A _ 1/2 ( K r)e lm<?i x 
— -~Jo 

xP-™ 1/2 ( C ose)e- im ^PrZ_ 1/2 (cosa), 







m— — oo 



(37) 
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FIG. 11: Analytic continuation to the imaginary axis of the 
angular momentum for the cone in the scalar case. 



where 9 and <f> are the angles of r in spherical coordinates, a and ip are the corresponding angles for fc, and is the 
angle between r and k. In the last equation, we have used the identity [1H| 



P„(cos0) = J2 (-l) m e im( ^ ' ) P I 7 m (cos6')P™(cos6i') , (6 + 9' < tt) 



(38) 



m— — oo 



where cos0 = cos#cos6>' + sin 9 sin 0' cos (cj) — cf>'). Equation (37) is only valid for Kc(9 + a) < tt/2, as can be seen 
from the asymptotic behavior of Legendre functions of a large degree [49], [50] , 



P l/ (cos( 



sin {(v + 1/2)6 + 7t/4) for large 



i/7T sm t 

By similar techniques, we can obtain the Green's function [49l [BTj . 

K 



Go{r,8,<l>;r',9 , ,<j)\iK) = — / cZA Atanh(A7r)fc lA _ 1/2 (Kr)/c iA _ 1/2 (fir')P A _i/ 2 (- cos 6), 



(39) 



(40) 



which in its bilinear form becomes 



+ 00 -oo 

Go(r,6,<t>;r',d',<f>',iK) =— y dA Atanh(A7r)fc lA _ 1/2 (Kr)fc lA _ 1/2 (Kr')x 



x e 



m— — oo ^ 

-(^^)p-™ /2(cose<) pm_ i/2( _ cose>)) 



where #< (0>) is the smaller (larger) of and 0' and we used the identity of Eq. (38) in the form [49 

oo 

P,(-cos0) = e 4m( ^*' ) P 1 7 m (cos6»)P™(-cos6'') 



> 



The Green's function can be also obtained by using the analog of Eq. ( 12 ) 

1 f°° 

dX Asinh(A7r)fc a _i/ 2 (r)fc iA _ 1/2 (r') = S(r-r'). 



(41) 



(42) 



(43) 



From the Green's function, we identify the scattering wavefunctions and the corresponding normalization factors, 

<&M,0) = k iX _ 1/2 ( K r)e im *Pr™ (cos8), 



C K \ m = — Atanh(A7r) 

47T 



(44) 



We note that Pj A ™i/ 2 ( C0S 9) is regular 5 everywhere except at 9 = tt (where it diverges) while P i ' A _ 1 / 2 (— cos 9) is regular 
at 9 = 7T but not 9 = 0, as we would expect for regular and outgoing wavefunctions, respectively. 



5 The superscript — m of the regular wavefunctions is chosen to absorb a complicated normalization factor. The distinction between the 
regular and outgoing wavefunctions arises from the arguments of the Legendre functions, not their order. 
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The T-matrices in the above basis can be easily found by matching boundary conditions. For a scalar field obeying 
Dirichlet and Neumann boundary conditions (labeled by D and N respectively) on a cone of half opening-angle 60, 
we have 

P- X m 1/2 (COS0 O ) 
° Xm ^- 1/2 (-cos0o)' 

TjVAm = "^^ f-COS0n) ' ( } 

Now we consider the cone opposite to an infinite plate obeying the same boundary conditions as the cone. We 



initially assume that the axis of the cone is perpendicular to the plate. From Eq. (37), the conversion matrix takes 
the form 

DxmMi = Atanh(A7r)(-l) m e- m ^^-i/ 2 (cosa) , (46) 

where fen is the wave vector parallel to the plate. We note that Re(a) = 0, so this expression is valid for 9 < ir/2. 
The T-matrix of the cone in the plane-wave basis is then 

/•oo +00 />oo poo +00 r-i 

T HK = \ ^ E / dX ' E ^Uta^. (47) 

7T7 — — (Yl U U 171 = — OO 



In the limit that 6 — > 7r/2, we can use the identity 

roc 

/ dX Atanh(A7r)^ 1/2 (coshi)F-™ /2 (coshu) = (-l) m 5(cosht - coshu) , (48) 
Jo 

to show that the T-matrix reduces to =F(27r) 2 <5(fc|| — fey), where the upper (lower) sign corresponds to Dirichlet 
(Neumann) boundary conditions. 
Then, for the Dirichlet case we have 

AC,w = - ~ A tanMA.) T <f X = m ± ^ ^ g A? M (2nd) , (49) 
2n 1 (iA + m + 1/2) "^-1/2^ cosOq) 



where 



2ir pioo 

A% mtX , m ,{x) = I di> I dcosa e'(---'^ i ,- ! n i/2(cosa)/ ,--' i/2 ( cosa)e -xcosa 



Ja=0 

oo 



/oo 
^A-i/2(0^A'-i/ 2 (0e- x « 



27r<5 m , ro , / dei^_ 1/2 (0^-i/ 2 (0e- xe . (50) 



We can then obtain the Casimir energy from 

he f°° 

£d = ^J o d,KtT]n(l Xm ,\'m' -■A/'Am.A'm')' ( 51 ) 
where the trace is defined as tr f\ m x> m ' — jfn° <iA 5Z^ = _ 0O fxm.Xm- The Neumann TV-matrix can be readily found by 



making the following substitution in Eq. ( 49 1 



P-_ 1/2 (cos0 o ) 4^T-i/ 2 (cos0o) 

(52J 



F-_ 1/2 (-cos0 o ) a|^A l -i/ 2 (-co^o) 

Next, we find an exact expression in the limit of small opening angle. Since the T-matrix vanishes in this limit 
for both the Dirichlet and Neumann case, we expand the logarithm to first order in the multiple reflections. For the 
Dirichlet case, we only need to keep m = term in the limit of small angle, since 

PAcOSdo) n 1 -.003), (53) 



P„(-cos0 o ) 21n(^) sin(i/7r) 
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and the energy becomes 

{^A-i/ 2 (cosa)P lA _ 1/2 (cosa)} + 0(0§) + • • • 

-- ^s\^m + ^ + -- (54) 

where the dots indicate corrections from higher reflections. 
For the Neumann case, in the limit of small angle, 



J-P„ m (cos0 o ) C rr ^inrnl^"^- t \ if W ^ 

ovq v v u/ J r(m)r(m+l)r(l+i>'— m) sin(i/-rr) V 2 / ) / 

^]_I^+1)(^)2 ifm = 

so we need to consider m = — 1, 0, 1. The Casimir Energy is then given by 

SfT = ^ f *« f d (cos a) e— r# fdAA^Mx 
W 2tt 8 Jo A ; Jo Jo cosh(ATr) 



(55) 



{-^A-i/ 2 (cos a )P iA _ 1/2 (cosa)(A 2 + 1/4) + 2^ A _ 1/2 (cosa)if A _ 1/2 (cosa)} + 0(flg) +••■ 

'"' ^o 2 +0(0 + ---, (56) 



24nd 

where in the last line we have used integral identities given in Ref. (51] . Again, the dots represent higher reflections. 

B. Electromagnetic field 

Because the cone does not have translation symmetry, the results for electromagnetism can no longer be obtained as 
a simple combination of two scalar problems. Using the same techniques as in the scalar case, however, we can obtain 
the electromagnetic Green's function for the cone in a useful form by analytic continuation of the usual representation 
in terms of spherical partial waves. As far as we know, this Green's function has not been obtained previously in the 
literature, so we derive it in detail in Section |Vf} We find that this case contains an additional subtlety, connected to 



the absence of the I = mode in electromagnetism (see Fig. 12). The Green's function takes the form 



Go = - £ £ f dX Atanh(A 7 r)^— (M-l 1/2im M^ 1/Vn - N° A ^ /2jm N^ 1/2 J 

oo 

E r(|m|)r(|m| + l)R™m R-o^* , (57) 

771 — — 00,771^0 

where we have defined the outgoing and regular wave functions 

M° A u i 1/2 , m = V x fc lA _ 1/2 ( K r)P™_ 1/2 (-cos0)e^r, 
M 7-i/2, n = V x k iX ^ /2 (Kr)P-^ /2 (cos0)e^r, 

N° A u l 1/2 , m = x V x fc jA _ 1/3 (Kr)i^_ 1/a (-co8d)e^r ) 

ri 

Rg£ = k (Kr) r x VP ~ |tn| (-cosfl)e im * , 

R^^oHrxVP Hra| (cos^ . (58) 

In this decomposition we have obtained the usual magnetic (transverse electric) modes M and electric (transverse 
magnetic) modes N, but we also have an additional set of discrete modes i?, arising from the pole at the origin in the 
contour integral. 
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FIG. 12: Analytic continuation to imaginary axis of the an- 
gular momentum for the EM case. 



From the Green's function, we can also read off the normalization coefficients 

M N K 1 

C\ m — ~C Xm = — — ^ 2 — ^ Atanh Air, 
C« = -^r(|m|)r(|m| + l). 
For a perfect reflector, the T-matrix is diagonal in the space of the M, N, and R modes, 6 

^A-l/2( C ° S ^) 
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R m 
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tan 
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2|m| 



We will also need the conversion matrices between the cone and plane wave bases, 
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E Oo 

m— — oo,m^0 
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(59) 



(60) 



(61) 



which can be obtained by analytic continuation of the expansion of a plane wave in vector spherical harmonics, 
yielding 

Mfey _ Atanh(A7r) 



^MXm,— A 2 + 1/4 ^iX-l/2\ LubLla )i 



L/4 

Mfey _ Atanh(ATr) 1 ^ (cmhn ) 

D " -~ X 2 + 1/4 ^sinh^ 6 ^i/ 2 (cosha), 



JVAm " 

Affey 
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ATfe„ 
MXm 



-|m|r(|m|)r(|m| + 1)^— e"^ P n Hm| (cosha), 

Slllll o 



-D 



Mfey 
AT Am' 



ATfe,, _ Mfc„ 
^NXm ~ MAro' 
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(62) 



6 In Fig. 1, M, JV and il modes are referred to as magnetic (M), electric (E) and the ghost (Gh) fields respectively. 
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where a — ia and ip = Zky are the angles of the plane wave in spherical coordinates, i.e., sinha = |k|| \ /k. 

We can now use these results to recast the T-matrix for the cone in the plane- wave basis appropriate for scattering 
from a plate. Symbolically and schematically, this transformation can be written as 

f°° p 1 

Tp, P , = £ / d\J2 CpD^ q j^-TqDq , (63) 

m Ja Q °<3 

where P, P' £ {M, N} are the wave functions in the plane wave basis whereas Q £ {M, N, R} arc defined in the 
conical basis. The A/"-matrix is then 

M Pk „ P% = JL^^^ T Pkll , P , k ,, (64) 



where 



P _ J -1 if P = M (Magnetic mode) 
[ 1 if P = N (Electric mode). 



We first consider a cone of small opening angle. In this limit, it is sufficient to consider m = of the mode N, in 

N 
N- 



which case only one component of the conversion matrix, D^, is nonzero. The T-matrix then simplifies to 



T N k^Nk' l} = y^W\J dX 1/4 cosh Att SinhQ p 'i\-i/2( cosha ) s ' mha ' P'ix-i/2( cosha ') + ( e o)- ( 65 ) 

To leading order, it suffices to consider only the first term of the multiple scattering expansion. Then to leading order 
as Oq the electromagnetic Casimir energy of a cone is given by 

= -2nj Q dK \ x dC ° Sha J d ^ j2^ e ~ 2KdC ° ShaTNk ^ + " 

he 1 1 f°° dcosha f°° tanhA7r 1 , .^,„ 2 > 

= -2^^^/ ^ho-io ' A cosh A7r + 1/4 ^-V 2 (^")^- 1/2 (cosh a) + P(flg) + • • ■ 

= -7 h ri,V °^'"- (66) 

a 167r | In ^ | 

where the dots represent corrections from higher reflections. This result has the same dependence on d and #o as the 
one for the Dirichlet cone at small angle 9 [see Eq. (54)], though the electromagnetic result is smaller by about 40%. 

For arbitrary opening angle, we compute the Casimir energy by keeping the first two terms in the multiple reflection 
expansion. The results are accurate within one percent and are not modified by higher order terms at the level of 
accuracy shown in Fig. 3. 



C. Tilted cone 



It is straightforward to allow the cone to tilt. Here we present the Casimir energy of a tilted needle — a cone in 
the limit of vanishing opening angle — opposite a plate. Since the T-matrix vanishes for small opening angle, we 
again keep only the first term in the multiple reflection expansion. To introduce a tilt, it suffices to consider the plane 



wave in a rotated basis. For the Dirichlet case, the conversion matrix is given by Eq. (46), where a and ip gi ye the 
spherical polar and azimuthal angles of fc, respectively, with respect to n, the symmetry axis of the cone. These are 
given in terms of spherical coordinates by 

cos a = cos a cos /? + sin a cos ip' sin j3, 
tan-0' 

tan^= -— — - (67) 

cos p — cot a' esc yr sin p 

where j3 is the tilt of the cone, and a' and tp' are the angular coordinates of k with respect to the x' , y' , z' axes, which 
form a convenient basis for the plate, see Fig. [13] These angles are defined according to 



In 2 + kl 

cos a' 



K 

ip' = Zk». (68) 
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FIG. 13: The configuration of a tilted needle opposite an infinite plane. 



Since a' is imaginary, we define a' — ia', so that 



cos a = cosh a' cos f3 + i sinh a' sin ip sin /3 
tan-0' 



tan ip 



cos /3 + i coth a' esc ip' sin /3 



(69) 



For the Dirichlet cone in the limit of vanishing opening angle — the Dirichlet needle — tilted at an angle j3 from the 
normal to the plate, the Casimir energy becomes, 



c needle 
£ n 



he 

2 



•c -1 f°° f°° 
— -r-, — 77. — 7777 I d« / d(cosha')e 
7r41n(V2)i h 



2nd cosh a.' 



dip' 



'. . tanh(A7r) 

COSh(A7Tj 



x ^a- 1/2 (cos a)P lA _ 1/2 (cos a*) + 0(6%) + 



he 



l 



l 



167rd|ln(6» /2)| cos^ 



O{0l) 



(70) 



Again the dots represent corrections from higher reflections. 

For electromagnetism, similar changes should be made in the conversion matrix to obtain the energy for a tilted 
needle, 



cnccdlc 
<~-EM 



he 



-1 



dn / d(cosha )e 



f\ „ — 2K,d cosh a 



27r41n(6» /2) 
x Pl x _ 1/2 (cma)Pl x _ 1/2 (cma*) + O(6 2 ) + ■ 

' 0(6l) + - 



dip' 



dX A 



tanh An 1 
cosh Xtt A 2 + 1/4' 



he 



1 



9(13)- 



I6ird |ln(6» /2)r^ y cos^ 



(71) 



where g(/3) is given by 



2lT j i f poo 



dip 

2^ 



d cosh a' 



cosh 2 i 



1 — cos a 



1 + cos a 



2tt j i l poo 



dlP 



da 



, tanh a' / (cosh a' cos pJ — 1 ) 2 + (sinh a' sin ip' sin / 



l 2ir J cosh a' y (cosh a' cos/? + l) 2 + (sinh a' sin ip' sin /3) 2 
The function g((3) goes monotonically from g(0) = ln4 — 1 to g(w/2) — 1 (see Fig. 3). 



(72) 



V. FINITE TEMPERATURE 



At finite temperatures the integral over the imaginary wave number k is replaced by a sum over Matsubara wavenum- 
bers, k„ = 2w i BT n. To gauge the importance of finite temperature effects, we report studies of two cases: 1) A 
perfectly conducting knife edge (a wedge in the limit of zero opening angle) opposite a conducting plate; and 2) A 
perfectly conducting needle (a cone in the limit of zero opening angle) either aligned vertically above a conducting 
plate or inclined at an angle p3 from the vertical. 
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A. Knife edge 



We consider a knife edge with an arbitrary tilt cf)Q. It would be straightforward to generalize the following results to 
nonzero opening angle. In Section [lll| we found the Casimir energy of this system by reducing the problem to a two 
dimensional one. To include the effects of finite temperature, the original frequency k must be restored and summed 
over Matsubara frequencies. We consider a single reflection (back and forth) between the two objects. As described 



in Section III this term gives a good approximation to the exact result. The free energy is given by 



F=-k B T 



2tt 



(73) 



i>0 



where T is the temperature and the primed sum indicates that the n = term is counted with a weight of 1/2. 
The dots represent higher terms in the multiple reflection expansion. For scalars subject to Dirichlet and Neumann 
boundary conditions, tr Af becomes 



knife r- -i 



1 



± 



\J K n + k z 



1 



V< + kl + kl \ y/nl + kl + k% cos^o 



(74) 



so that for electromagnetism the first term cancels between Dirichlet and Neumann modes and we get an overall factor 
of 2 (summing over both modes). This yields the EM Casimir free energy to first order in the multiple scattering 
expansion, 



knife 

edge 
EM 



hcL, 



1 



1 d cQth d 
16-7T 2 cos 4>q d 2 At At 



where At 



He 



2nk B T 



Expanding the force at low temperature, we find 



knife 
podge 
r EM 



hcL, 



1 1 



8ir 2 cos 0o d 3 



1 

45 V A T 




(75) 



(76) 



Here we have found the contribution from the first reflection which shows that the force depends only weakly on 
T at low temperature. The dots represent the contribution from higher reflections, but as argued before, the latter 
is significantly smaller. It is interesting to note that the last equations are quite similar to the result from two 
parallel plates. To the first order in multiple reflections, the free energy of two parallel perfectly reflecting plates at 
temperature T is 



J", 



parallel 
plates 
EM 



, ^L X L Z d ( \ , d 
k sT—^ — i-jcoth — 
8tt ad \d \ T 



(77) 



When e = ir/2 — <fio is small, the result for the wedge agrees with the result for parallel plates, if e is identified correctly 
in terms of temperature and other parameters of the problem. 

In contrast, our numerical results follwing from Eq. (74 1 indicate that thermal corrections for a scalar field obeying 



Dirichlet or Neumann boundary conditions on a knife edge are proportional to temperature with a power close to 
3. This points to the complex interplay of temperature dependence and geometry, which was previously noted in 
Ref. I321. 



B. Needle 



While temperature corrections can be computed for cones of arbitrary opening angle, we can obtain a closed-form 
analytical formula for a needle. Hence we focus on the latter case below. Our calculation exploits the fact that the 
T-matrix tends to zero for a sharp needle and thus the leading term in the multiple reflection expansion is sufficient. 
We will see that at any nonzero temperature, the free energy diverges in the infrared. This divergence doesn't depend 
on the objects' separation, however, so the force remains finite. 

To first order in multiple reflections, the free energy is 



T = k B T^2' trAT- 



(78) 



K„>0 
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For Dirichlet boundary condition, this gives 



7 rp ^ />00 -i 

J -needle = _ «B ^ / da; eX p(-2^) + 0(^ 2 ) + • • • . (79) 

From the last equation, it is clear that the energy becomes infinite for any finite T due to the n = term of the sum. 
However, if we take the derivative with respect to d to get the force, we obtain 

7 rp OO ' />00 

F "° dIC = ~ 4|ln1>o/2| ^ 2Kn I dX ex P(" 2K « da; ) + °(*o) + ' ' ' ■ (80) 

Note that the integral must be computed first; only then can the sum be performed. If we separate the summation 
into the first term (n — 0) and a sum over all n > 1, however, the latter sum commutes with the integral. We then 
obtain 



This is equivalent to Eq. (54) for T —> 0, while it becomes linear in T for T —> oo. An interesting feature is the 
leading thermal correction for small temperatures. By expanding the last expression, we find 

he 1 



TTinccdlc 



167r|ln0o/2| d 2 



lfd\ 2 ,//d^ 4 



3 \ At / IV 



O(0 2 ) + ■■■. (82) 



This shows that thermal corrections to the force have a significantly larger effect than in all other known examples, 
most notably two parallel plates 7 . 

For electromagnetism, we find the result 



Again this is equivalent to Eq. (66) for T — > and it scales linearly with T at large temperature. Expanding this 
equation at low temperature, we obtain the leading thermal correction, 

*ST = - 4 H2d/X T ) + 0M0 2- + • • • ) + OiSl) + • • • , (84) 

where the coefficient of the third term equals |(— 1 — 27 — In 4 + 24 In A), in which 7 is Euler's constant and A the 
Glaisher-Kinkelin constant. Interestingly, the first thermal correction is again proportional to a (second) power of T 
but now also multiplied by the logarithm of d/Xr- 

We can extend these results for a needle tilted by an angle f3. The Dirichlet free energy is just multiplied by 1/ cos/3, 

tilted i t 1 

F - CdlC - -^lnt72M C ° th ^ W + °W + ^ 

For electromagnetism, the /3-dependence changes to 



^needle. 9 {fi,d/\ T ) 2 

F ™ -^167r|lnV2M 2 cos {3 + ° {do) + ' " " ' (86) 



where g((3,d/\T) is defined as 



±,fd_ \ 2 p / (xcos/?-l) 2 + (:r 2 -l)sin 2 /3sinV 1 

5^,0/ at; At + ^ At J y Q 2?r ^ V (x cos /? + l) 2 + (.x 2 — 1) sin 2 f3 sin 2 sinh 2 (o?/ \t)x ' [ ' 



7 Relatively large thermal corrections, O (T 3 ), have been reported for Dirichlet half-plate geometry in Ref. 1321 . 
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From this equation, the thermal correction at low temperature changes from 0(T 2 InT) for a vertical (EM) needle to 
0(T 2 ) when the needle becomes almost parallel to the plane. Eq. (84) and (87) predict large temperature corrections 
to the Casimir force at room temperature, as shown in Fig. 3. 

The striking dependence on geometry of both the analytic form and the size of finite temperature corrections was 
suggested in Ref. [S^ and studied for a scalar field obeying Dirichlet boundary conditions in Ref. [35]. Here we 
found that thermal corrections are relatively small for a wedge/knife-edge (along with parallel plates) while they are 
significantly larger for a needle. 



VI. IMPERFECT CONDUCTIVITY 

Throughout this paper, we assumed that the objects are perfect reflectors for the EM field. In this section, we 
estimate to what extent this is justified for realistic conductors. Our analysis is based on the fact that the contribution 
to the Casimir force between two objects is suppressed at wave numbers greater than their inverse separation. Small 
frequencies contribute most, just where the assumption of perfect conductivity is best. Although the distance scale 
for the wedge-plate and cone-plate geometries is d — the tip (or edge) to plate separation — the typical separations 
that contribute to the Casimir force for these open geometries is much greater than for parallel plates. 

To investigate this further, we write the Casimir energy as an integral over (imaginary) wave number, 



£ = 



dn p(n, d) , 



and compare the Casimir energy density, p(«, d), for the cone-plate geometry with the parallel plate case. For a 
perfectly reflecting needle opposite an infinite plate, this density function can be read off the first line of Eq. (66), 



p(K, d) ~ C I dp 



P 



1 1 



p + 1 p 



(89) 



where C is a constant which depends on the opening angle. Plots of the energy density are given in Fig. |14| It is 
clear that the energy density of a needle-plate system is falling off rapidly as opposed to parallel plates where it falls 
off rather slowly. 




2nd 



FIG. 14: Energy density p(K,d). The curves are not normalized and should be scaled with appropriate parameters of cone or 
plate geometry. The dashed line represents the energy density for two perfect metal plates, the solid line the energy density 
for a perfect metal cone and a plate. It is obvious that small frequencies are more important than large ones for the cone-plate 
as compared to the two-plate case. 



To be more quantitative, we examine the Casimir energy as a function of the upper limit on the wave number 
integral, 



£ (k, d) 



dn' p(k , d) , 



(90) 



and from this define a maximum wave number, K maX j above which the contribution to the Casimir force is negligible. 
The plots of £ (k, d) normalized by £(oo, d) are given in Fig. 15 for the same two systems. The cutoff, K max , is defined 
by 



£(K max ,d) = 95% £{oo,d), 



(91) 
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which gives an upper bound on the frequency at the price of 5% inaccuracy. The cutoff can be read off the plots which 
give Kmax ~ 0.5/d for a needle-plate system and K max ~ 2/d for two parallel plates. The cutoff for the cone-plate 
system is considerably smaller, as anticipated. 

£(k, d)/£(oo, d) 




FIG. 15: Energy £ (re, d), normalized by £(oo,d) (where, as usual, all frequencies are integrated over), is plotted as a function 
of red. The dashed line pertains to two perfect metal plates, the solid line to a perfect metal cone and plate. The energy in 
the cone-plate configuration is dominated by contributions from small frequencies when compared to the two-plate case. 

So far, we considered a perfect reflector and required that 95% of the energy is contributed by frequencies smaller 
than K max . Next we address the question whether this regime can be approximated by perfect reflectivity. A finite 
dielectric function only affects the T-matrix. Here, we derive an integral equation which determines the T-matrix and 
allows us to study the deviation from perfect reflectivity in a systematic way. 

In the limit of a sharp cone, it suffices to consider only the m = components of the field. There is no R mode in 
this sector. Hence, we need to solve the matrix equation for the T-matrix only for electromagnetic M and N modes 
at m — 0. A further simplification occurs for m = since the polarizations decouple. Finally to lowest order in 
the opening angle, the M mode can be neglected, as it can for the case of perfect reflectors. In summary, we must 
compute the scattering of the N mode only at m = 0. 

By imposing the continuity conditions for electromagnetic field, we find for the T-matrix elements with m = the 
integral equation 

I 

d(3 T afs d^ (P i p_ 1 / 2 (-cos9o)de o P il - 1 /2(cos0 o ) - -^ P^_i/ 2 (- cos0 o )-Pi 7 -i/2(cos0o)) 

=d ai (P M _ 1/2 (cos6>o)d eo P, i7 _ 1/2 (cos6»o) - ^d 6o P ia _ 1/2 (cos 9 Q )P n _ 1/2 (cos O )) , (92) 
where e is the permittivity and 

/■OO 

dfa= dx k. l p_ 1/2 (x)k n _ 1/2 (nx) 1 (93) 
Jo 

with n = \J~k. If e is large, we can systematically organize the perturbation theory as 

T = T° + T 1 + ■ ■ ■ , (94) 

where each order is smaller by a factor of 1/e and T° is the perfect reflector result. Using this ansatz in the matrix 
equation and exploiting the smallness of #0j we obtain 

T «f> = - e (iL) W^klW col?J dl da -<^Tyl ( °^' (95) 

Although we will not compute this integral explicitly, it can be compared parametrically with the zeroth order result, 
T° ~ 1/| ln#o/2|. The first order solution is smaller than the perfect reflector result parametrically by a factor of 

(96) 



e(icK) <92|ln(9 /2| 



Of course, for extremely small angles, perfect reflectivity is incorrect. However, if the permittivity e is large enough, 
the correction can be small, even at small angles. Note that, at worst, the frequency in the argument of e should be 
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replaced by K max , which is given by 0.5 jd for a sharp cone. Using the plasma model 

*->- 1+ S&- (97) 

for gold (A p = 137nm), and choosing a relatively small angle 9q = 0.1 at a fairly small separation distance of 
d = 200nm, we get 5 5%. The correction to the cone's T-matrix implies the same correction to the energy (in the 
first reflection, the energy is just proportional to the T-matrix). Therefore, the assumption of perfect reflectivity 
seems to be quite good even for small angles and short distances. Here, we studied the finite conductivity corrections 
only for the cone to estimate its importance. The same effects for the infinite plate can be found exactly within 
this formalism since the T-matrix of a dielectric plate remains diagonal. By employing the same dielectric model at 
the same separation distance (d — 200wn), we can easily see that finite conductivity correction for an infinite plate 
(facing a perfect conducting needle) is just a few percent. 



Appendix: The electromagnetic Green's function 



In this section we derive the electromagnetic Green's function in coordinates appropriate to the cone problem. 
Similar, though more tedious, algebra will reproduce Eq. (62 1, which gives the conversion matrices. 
We start from the dyadic Green's function in terms of spherical vector wave functions, 

oo oo 

G Q (r;r', ik) = K^T ^ V x i t ( K r <)¥{"■ (r)r ® V x fy(KT>)Y; m * (r')r' 

I — I 771 — — OO 

- -V x V x iAnr^Y," 1 (r)r® -V' x V' x kAnryW," 1 * (y')y' 

K K 

and write this expression in a form which is more appropriate for our purpose, as 

oo l „, 

G (r;r', iK )=£ £ „ f ^ *+ (Mjfffr) ® M£T (r) - Nff fr) ® N^* (,)) , 



2=0 m= — l 



1(1+1) 47T 



with 



M™ g = V x i, («r)if* (cos 0)e im *r, 
N ta? = -V x V x i l (nr)P{ n (cos9)e lm ' l 'r. 

The outgoing wave functions are defined similarly, with kg{nr) substituted for i^{kt). Then, 

M ?m <^ M L S * " N °m ® N [m* = D [H^P^ (cos ff)^' (-\) m %i (/sr)P,- m (cos 9)e 
where the dyadic differential operator D is defined as 



V x r (g) V' x r' - V x V x r ® -V' x V' x r'. 

K K 



The summation over m in Eq. (99 1 can be extended to (—00,00) because 

PI™ (cos 8') PC m (cos 9) = 0, for I < \m\ and I, m integers. 
By using P™(cos0) = (— l)' _m P™(— cos 0), we can write the Green's function as 



00 00 

l— l m— — 00 



1 2Z + 1 

Z(Z + 1) 4tt 



(-1)' k l (Kr')P l m (-cosO')e im<1 '' i,(«r)Pf m ( 
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(98) 



(99) 



(100) 



(101) 



(102) 



(103) 



(104) 



In order to turn the sum into a contour integral, we should divide it by sini/7r to obtain poles at the integers v = I. 
However, unlike the scalar case, we now have to exclude I = from the sum. Moreover, the factor of l/l in front of 
the sum defines a pole at the origin. Therefore, analytic continuation of the integral yields two different contributions: 
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an integral from —1/2 — ioo to — 1/2 + ioo denoted by C , and a small circle at the origin denoted by Cq (see Fig. 12 1. 
From the asymptotic behavior of the Legendre functions given in Eq. (39), we see that the large half-circle at infinity 
can be neglected for & > 9. So we have 



where 



"0 



m— — oo ' 



with the vector wavefunctions defined in Eq. ( 58 1 , and the residue of the pole at the origin gives 

dv- 



1 2V + 1 7T 



L^Co 



z^(y + 1) 4-7T sin7r^ 



k v {Kr')P™{- cose 1 ) i„{Kr)p- m {cos9)e- im ^-*' ) 



(105) 



(106) 



(107) 



Because of the double pole in the integrand, we should take the first derivative of its coefficient (with respect to v). 
This is further simplified by noting that -fo"(~ cos 9')P^~ m (cos 9) = (except when m = 0, which does not contribute 
to the Green function since it vanishes upon applying the D operator). So we must compute 



d u (Pl n (-cos9')p- m (cos9)) 



A bit of algebra gives 



"■Co 



Y fc o (^')jo(^)r(|m|)r(|m| + l)P ( 7 |m| (-cos0') P o - |m| (cos0)e 

47T — ' 
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(108) 



(109) 



This is not yet in the form that we need, because the radial coordinates r and r' should appear symmetrically. The 
symmetry becomes manifest by applying D. The double and quadruple curls in the D operator combine to yield 



G 



Co 



£ )T T(\m\)T(\m\ + 1)R°^ R\ 



4tt 



out -p rcg* 
0,m 



(110) 



m— — oo,m^0 



where the vector functions R are defined in Eq. ( 58 1 . We thus obtain the full Green's function, Eq. ( 57 ) . 
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